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1. INTRODUCTION 

The feedback linearization technique is one of modern tools, which allows synthesizing a control 
law for smooth, continuous, non-linear systems [1-5]. The basic idea is to transform a nonlinear system into a 
linear one via nonlinear change of coordinates only (state linearization problem) or by nonlinear feedback 
and change of coordinates (feedback linearization problem) so that the linear control techniques can be 
applied [1]. Feedback linearization has been one of the most active research topics in recent years; 
linearization of affine systems [6], systems with nonsmooth nonlinearities [7], exact linearization [8]. Various 
methods exist to find feedback linearizable form for single input [9], and multi-input nonlinear systems [10-15]. 

Significant research effort has been devoted to the construction of approximate solutions to the 
problem of linearizing nonlinear systems by state or output feedback. Main reason for that is the limited 
applicability of the rigorous methods, and the complexity, sensitivity and design difficulties of the exact 
linearizing compensators, if any [16]. The analysis and the comparison of approximate linearization and 
exact linearization are proposed by [17]. If the original nonlinear system cannot be linearized exactly by state 
feedback, the method of approximate feedback linearization [16, 18-22] is used. This process consists of 
finding approximate output functions that satisfy the involutivity condition up to determined system order, 
[16, 19] adopted a two-step procedure to solve the approximate linearization; first a state transformation 
matrix is settled, so that the nonlinear system is transformed approximately into the controllable canonical 
form. Second, a standard nonlinear linearization method is used to transform the controllable canonical form 
into a stable linear system. 
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Kabanov [19] offered an approach that consists in representing the original nonlinear (single and 
multivariable) system into a state-dependent coefficient form and applying similarity transformations that 
allow getting the system to canonical form, that considerably simplifies the problem of feedback linearization 
afterwards. Such similarity transformations allow accomplishing linearization of system without considering 
the virtual output function y = h(x) [19].This work consists in applying a block companion similarity 
transformation to a class of nonlinear multivariable system; where the number of inputs divides exactly the 
number of states. It is known that if a linear-time-invariant MIMO system described by a state space equation 
has a number of states divisible by the number of inputs and it can be transformed to block controller form, 
we can design a state feedback controller using block pole placement technique by assigning a set of desired 
Block poles. These may be left or right block poles [23-26]. The eigenvalues and eigenvectors 
(eigenstructure) of the state matrix can determine system performance and robustness more directly and 
explicitly than other indicators. Hence their assignment should improve feedback system performance and 
robustness distinctly and effectively [25]. 

In the present work, firstly we have started by introducing feedback linearization used in 
multivariable systems. As a second section a problem of transforming a multivariable system into block 
controller form is formulated, it is then followed by the discussion of the feedback linearization using 
similarity block transformation for the class of MIMO system which is the heart of this work. In a forth 
section examples are used to show different steps and the effectiveness of the proposed method, and finally 
the paper is finished by a conclusion. 


2. PROBLEM FORMULATION 
Consider a class of multivariable nonlinear system 


x = f(x) + X1 gi) ui 


x= f(x) + G(w)u (1-a) 
y, =hi(x),i=1,..,m (1-b) 
or Y = H(x) 


where f, gand h are sufficiently smooth, m is the number of inputs and outputs,n is the number of states. 
The multivariable system (1), in the state-dependent coefficient form is given by: 


x(t) = A(x)x + B(x)u (2) 


where A(x) is ann X n and B(x) is ann X m matrices respectively; are continuous differentiable. 

Consider the system (1), the feedback linearization problem consists in finding the nonsingular 
coordinate transformation Z = Tp; (x) and the input control law u = ug, such that the system (1) can be 
transformed into linear block canonical form. Under condition of the existence of smooth function H(x), 
transformation Z = Tp; (x) is described by (1). 


H(x) 
Tno) =| FHC 6) 
LH (x) 
and linearizing control is given as, 
Up, = S*(x)(v — q(x)) (4) 


where q(x) = LH (x) and s(x) = Lg Ly *H (x)and v is the new control input;L-H (x) is the Lie derivative of 
the it” order from the function H(x) along the vector field f [1, 19]. 

Applying feedback linearizing control ug, (4), in respect of the transformation Tp,(x) (3), the 
system (1) has the following form, 


4=A,2 + B,v (5) 
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where (A,, Be) is in block canonical form. 


aa Om 

as ei ~ E [B= lo (6) 
| rn = 
Oe, 20s, m 


where 0,, and Im are m X m zero and identity matrices, respectively. 
Applicability of the normal form method for feedback linearization is ensured by the existence of 
diffeomorphism (3) and observance of the terms [1, 19], 


LGH (x) = LgLpH(x) = ... = Lely 7H(x) = Om 
s(x) = LeLF'H (x) # Om 


which is not always true. In such situations, as a rule, methods of approximate feedback linearization are used 
[16, 19, 27]. The block canonical similarity transformation used for feedback linearizing control design is: 


z = T(x)x, x = T1 (x)z,T (x) E R™” (7) 
3. FEEDBACK LINEARIZATION 

System (1) can be represented in the state dependent coefficient form [28]: 

x(t) = A(x)x + B(x)u (8) 
Let us assume that matrix A(x) is n X n and vector B(x) is, n X m is continuously differentiable. It is also 
supposed that the pair (A(x), B(x)) is a full block controllable, i.e. for the system (8) controllability matrix 
has full rank (rank of w(x) = J), for t = to. 


w(x) =[B(x) A(x)B(x) .. AUt(x)B(x)] 


where = = l, l (an integer) is the number of blocks. 


Applying transformation (7) to (8) we get (5) for which the following relations are true. 
z(t) =T(x)x + T(x)x 


2(t) = T(x)T1(x)z + T(x) [A(X)T “1 (x) z + T(x) B(x)u] 


Z(t) = A,(z)z + B.@)u(z) (9) 
A(Z) = T(x)A(x)T “1 (x) + POT) (10) 
where 
Om la a Om 
A-=| : a (11-a) 
: : : en 
—Ap(z) —-Ai(Z) = —Aj-1(Z) 
Om 
B, = T(x)B(x) = i (11-b) 
m 
Im 
where A; E R™™,i = 1,2,...,.1. Im and Omare m x m identity and null matrices respectively, | is an integer 


such that | = = where n is the order of the system and m is the number of inputs. 
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The right matrix fraction description (RMFD) [24-26] of the system can be formulated directly from 
(11-a) as, 


F(A) = Npe(A)Dg* (A) (12) 
where, 

DR(A) = Ima! + Aad 14... +4, (13) 

Ng(A) = IA! + CA’ 34... +C, (14) 


Ci E R™™ are constant matrices. The desired matrix polynomial Ap(A)if R,is a right solvent is given by, 
Ap(A) = R! + DiR! + + D, 
[Da Daa—1) -Dar | = —[RIR} ... Ri] Vg 


where Vpis the right block Vandermonde matrix[24, 25]. An alternative factorization of F(A) is the left 
matrix fraction description (LMFD) defined by the following equation; 


F(A) = D ANA) (15) 


where D,(A) is the left denominator (13) and N,(A) is the left numerator (14). The desired matrix 
polynomialA p (A) if L; is a left solvent is: 


Ap(A) = L + LtD, +++ D, 


the desired solvents are given as, 


Darl 
Da Im bi a IT HE 
Dau-1)| __|Im L2 = L51) [ib 
Da Tees Vilage Sao SE EWE 


Linearizing control for block companion system (11) can be defined in the form of the feedback; 
u(z) = Up, (Zz) + v(z) = K(z)z + G,(z)z (16) 


where matrix G,(z) can be calculated, for instance, by means of the method of the placement of the closed- 
loop system block poles. K (z) is computed by [19, 21]. 


K(z) = —[Ara@) AZ) -. Ao(z)] (17) 
In terms of initial system state variables the linearizing control (16) can be written as, 

u(x) = Up, (x) + v(x) = K(x)T (x) x + G.T (18) 
The procedure to find the coordinate of change is [24, 25]: 

Tx) = {Tix}, i = 1,...,1 (19-a) 
let 


Teor = [0m Om e Imjwo*(x) (19-b) 
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Ter 


Tc 


T(x) = (19-c) 


T, A171 


where Omand Im are null and identity matrices, respectively. The transformation (7), (19) exists, if and only 
if: 

- controllability matrix w(x) has full rank n at t > to 

- = = l,l is an integer 


If we suppose that the block canonical form has the form (6). 


On Om On Om 
Ac — : : š B: = 0. 
O Wia Ei 
Om Om Om c 
ie. K(z) = [0m a On] (20) 
The linearizing control is: 
u(z) = v(z) = G,(z)z = K.(z)z (21) 


In the next section, numerical examples are illustrated to show the effectiveness of the proposed method. 


4. ILLUSTRATIVE EXAMPLES 
4.1. Example 1 
Let consider the following nonlinear multivariable system, with a number of states is n = 3 and a 


number of inputs is m = 3 therefore a number of blocks is l = = 1. The control problem is to stabilize the 


system (22) by means of the feedback. 
Ky = x1 +x tu, 


Žž = —x1 + 2x3 + 2u, (22) 


X3 (x pg 1)x2 + U3 


where, 
x, +x Uy 
f(x) =|—x, + 2x, | and g(x) = | 
(x, — 1)x2 Us 


We write f(x) and g(x) in state-dependent coefficient form; A(x) and B(x), respectively. 





1 x 0 1 0 0 
A(x) = f- 0 =2|,B(x)= f 2 o| (23) 
x2 —1 0 0 0 1 
We compute the similarity transformation matrix. 
1 0 0 
w(x) = [B(x)] = f 2 o| (24) 
0 0 1 


The rank of w(x) is 3 (full rank), hence the system is controllable. Now we compute the coordinate of 
change T(x)we calculate first. 
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Ta = [Ip]w* e) 


then, 
1 0 0 
T@)=[a]=|]0 5 0 (25) 
0 0 
hence, 
A, (z) = T(x)A(x)T™ (x) 
1 2x, 0 
A(Z) = -ż 2u 24 (26) 
X% —2 0 
and K(z) = —A,(z), and we have 
1 0 0 
Bo = 769800) =[o 1 o| (27) 
0 0 1 


In this application, G,(z) is computed using desired block poles of the closed loop system. The 
eigenvalues of close-loop system matrix are 2, = —3, Az = —4 and A3 = —5. The characteristic matrix 
polynomial of the closed-loop system is forced to equal a desired matrix polynomial of the form. 


Ap(A) = Ima! + DA 1 +- + D; (28) 


The desired characteristic matrix polynomial Ap(A)may be constructed from a set of desired m x m 
block poles or matrix roots of Ap(A); These block poles can be constructed under different forms either an 
m Xm diagonal form matrix, an m X m controller canonical form matrix or an m X m observer canonical 
matrix. Using block poles in diagonal form gives a better response, robustness and smaller sensitivities [24], 
for this reason only the case of the diagonal form will be studied. Feedback gain of the closed-loop system 
G,is computed as: 


A-— BG. =AÁp (29) 


Ap is a desired closed-loop matrix in block diagonal form where each block is constructed from 
desired eigenvalues. The desired block pole constructed in diagonal form is given by, 


-3 0 0 
R=|0 —4 0 | (30) 
0 0 -5 





The corresponding 3 x 3 desired right denominator matrix polynomial of degree 1 is: 


3 0 0 
Dr(s) = Is + Djıwhere Dp, = —R?Vg' i.e., Dey = f 4 o| 
0 0 5 


Now we can find Ge 


3 2x, 0 
G.=|-"/, 4 -1 
X2 —2 5 


The linearizing control. 


u(z) = Up, (zZ) + v(z) = K(z)z+G,(z)z = [K(z) + G,(z)]z = K,(z)z 
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where 
-3 0 0 
Kc=|0 —4 0 | (31) 
0 0 —-5 





its norm ||K,||, = 5. In the terms of initial system state variables linearizing control (18) can be written as, 


u(x) = K(x)T(x)x + G(x) T (x) x = (K(x) + G,(x))T (x) x = K,(x)T (x) x 


hence, 
-3 0 0 
u(x) = | 0 -4 0 | (32) 
0 0 —-5 


The results of the nonlinear system (22) with the linearizing control (32), For modeling the following initial 
values were taken: x,(0) = 0, x (0) =1,x3(0) = —1. Using MATLAB software (M. file), the state 
variables responses are shown in this Figure 1. The results of the figure demonstrate that the system is stable. 



































-0.4 





-0.6 











-0.8 |- 

















o 0.5 1 1.5 2 
time (s) 


Figure 1. State variables result for the system (22) 


4.2. Example 2 
In the following nonlinear multivariable system, a number of states is n = 4 and a number of inputs 


is m = 2 therefore a number of blocks is l = - = 2. The control problem is to stabilize the system (33) by 
means of the feedback. 


ý = x +x + x2 

ž = —X1 — 2X3 + Uy, (33) 
Xz = X4 

Xq = (x1 — 1)x, + u, 


where, 
X +x + x3 0 
fe) =| “> 2%8 | and g(x) = | 
X4 0 
(x, — 1)x2 Ug 


We write f(x) and g(x) in state-dependent coefficient form; A(x) and B(x), respectively. 


1 1 x 0 0 0 
_|-1 0 -2 0 _|1 0 
x% —-1 0 0 0 1 
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We compute the block similarity transformation matrix. 


0 
1 
0 


ooo 
oor 


w(x) = [B(x) A(x)B (x)] = (35) 


0 
0 
1 
0 1 -1 0 


The rank of w(x) is 4 (full rank), hence the system is controllable. Now we compute the coordinate 
of change T(x), we compute first. 


0 0 1 0 


Ta =[9 9 o ">e 
then, 
10 0 0 
rospi tay 0 G6) 
0001 
hence, 
0 0 1 0 
a ETAT OE 4 > ea, 67 
x2+1 x3; -1 0 
and we have 
0 0 
Be(z) =TH)BO) =|) g (38) 
01 


In this application, G,(z) is computed using desired block poles of the closed loop system. The eigenvalues 
of close-loop system matrix are A, = —1, à, = —2 and A, = —3 and A, = —4. The desired block pole 
constructed in diagonal form is given by, 


_f-1 0 _f[-3 0 
R=[; ` | anar = [e Al (39) 
hence, Dy(s) = Is? + Djs + Dro 
The linearizing control. 
u(z) = up, (z) + v(z) = K(z)z+G,(z)z = [K(z) + G,(z)]z = K,(z)z 
where, 
f2 —2 5 X3 
ke=| 4 x3+8 -1 6 GD 
In the terms of initial system state variables the linearizing control is. 
_ [0.43 14 1.4x3;-2 x3 
HOS | x —1 008 06 (41) 


The results of the nonlinear system (33) with the linearizing control (41), for modeling the following 
initial values were taken: xı (0) = 0 , x (0) = 1,x3(0) = —1 and x,(0) = —0.5. Using MATLAB software 
(M. file), the state variables responses are shown in this Figure 2. The system is stable as it is indicated in the 
results of the figure. 
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Figure 2. State variables result for the system (33) 


5. CONCLUSION 

In this work, state feedback linearization of a class of nonlinear multivariable systems; the number 
of states divides exactly the number of inputs, by using block canonical similarity transformation is 
presented. First the original nonlinear system is represented in state-dependent coefficient form then a block 
companion similarity transformation is applied to the system. Linearizing control for block canonical system 
is defined in the form of the feedback; where the feedback gain is computed by forcing the characteristic 
matrix polynomial of the closed-loop system to equal a desired matrix polynomial constructed from a set of 
desired block poles or matrix roots in diagonal form. In this approach some advantages are taken: i) Because 
there exists an approximate linearization method when the exact feedback linearization is not applicable, ii) 
Because using block controllable form, the eigenvalues and eigenvectors (eigenstructure) of the state matrix, 
can determine system performance and robustness, iii) Using block poles (solvents) in diagonal form gives a 
better time response, robustness and smaller sensitivities of the system, iv) Such block similarity 
transformations allow accomplishing linearization of system without considering the virtual output function 
Y = H(x) or y; = h;(x). In this work, similarity block transformation is developed for class of MIMO 
system, which seems to be a successful approach. As further work, an enhancement to the present method is 
proposed, the similarity block transformation will be applied to a general form of MIMO systems; for the 


ny. . 
case where l (or —) is not an integer. 
m 
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